The Onset of Quark-Hadron Duality in Pion Electroproduction 
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A large data set of charged-pion (tt ) electroproduction from both hydrogen and deuterium tar- 
gets has been obtained spanning the low-energy residual-mass region. These data conclusively show 
the onset of the quark-hadron duality phenomenon, as predicted for high-energy hadron electropro- 
duction. We construct several ratios from these data to exhibit the relation of this phenomenon 
to the high-energy factorization ansatz of electron-quark scattering and subsequent quark — » pion 
production mechanisms. 

PACS numbers: 12.40.Nn, 13.87.Fh, 12. 39. St, 13.60.Lc 



At high energies, the property of Quantum Chromodynamics (QCD) known as asymptotic freedom allows for an 
efficient description in terms of quarks and gluons — or partons, weakly interacting at short distances. In contrast, 
at low energies the effects of confinement impose a more efficient description in terms of collective degrees of freedom, 
the physical mesons and baryons — or hadrons. 

Despite this apparent dichotomy, in nature there exist instances where low-energy hadronic phenomena, averaged 
over appropriate energy intervals closely resemble those at asymptotically high energies, calculated in terms of 
quark-gluon degrees of freedom. This is referred to as quark-hadron duality, and reflects the relationship between the 
strong and weak interaction limits of QCD — confinement and asymptotic freedom. 

The observation of this phenomenon in fact preceded QCD by a decade or so, with remarkable similarity found 
between the low-energy cross sections and high-energy behavior in hadronic reactions, with the former on average 
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appearing to mimic features of the latter. At that time, this was explained with the development of Finite Energy 
Sum Rules, relating dispersion integrals over resonance amplitudes at low energies to Regge parameters describing 
the high-energy scattering Q. The equivalence, on average, of hadron production in electron-positron annihilation 
and the underlying quark-antiquark production mechanism was later similarly understood Q . 

It was natural, therefore, that this same framework was used to interpret the early observation of quark-hadron 
duality in inclusive electron-nucleon scattering. Bloom & Oilman found that by averaging the proton Fi structure 
function data over an appropriate energy range, the resulting structure function in the resonance region closely 
resembled the scaling function which described the high-energy scattering of electrons from point-like partons [4J. 
Recently, the phenomenon has been revisited with unprecedented precision, and was found to work quantitatively far 
better, and far more locally, than could have been expected 

Although postulated to be a general property of QCD, the dynamical origin of quark-hadron duality remains poorly 
understood. It should manifest itself in a wide variety of processes and observables. In this Letter, we generalize the 
duality concept to the unexplored region of ( "semi- inclusive" ) pion electroproduction eN — > eir^X, in which a 

charged pion is detected in coincidence with a scattered electron. The missing mass of the residual system X, M x , 
is in the nucleon resonance region (defined here as M 2 < 4 GeV 2 ) for the remainder of this Letter, and we will show 
the dual behavior of this region with a high-energy parton description. 

At high energies, perturbative QCD predicts factorization between the virtual photon-quark interaction and the 
subsequent quark hadronization, 

dn^dE^dzdp'^dct, dN _ bp 2 1 + Acos{(j)) + Bcos(2(j)) 

da = ~dz T 2tt ' ' ' 

dn e dE c , az Z7r 

^-^e^q^Q 2 ) D q ^(z,Q 2 ), (2) 
q 

where the fragmentation function D q ^ 7T (z,Q 2 ) gives the probability for a quark to evolve into a pion ir detected 
with a fraction z of the quark (or virtual photon) energy, z = E w /v. The parton distribution functions q(x,Q 2 ) are 
the usual functions depending on the Bjorken variable x and Q 2 . The transverse momentum p?, z and the angle 4> 
reflect the extra kinematical degree of freedom associated with the pion momentum. Both the parton distribution 
functions and the fragmentation functions depend on Q 2 through logarithmic Q 2 evolution. Their dependence on 
Pt is removed in a Gaussian approximation, reflected in the noted exponential px dependence, with b the average 
transverse momentum of the struck quark. In the (very) high energy limit, the factors A and B become zero. At lower 
energies, these "factors" reflect the longitudinal-transverse and transverse-transverse interference structure functions 
of the general pion electroproduction framework [9(, and can, e.g., vary with z and Q 2 . Note that a consequence of 
this factorization ansatz is that the fragmentation function is independent of x, and the parton distribution function 
is independent of z. 

At lower energies, where hadronic phenomena dominate, it is certainly not obvious that the pion electroproduction 
process factorizes in the same manner as in Eq. @. However, it has been argued that at relatively low, yet sufficiently 
high energies, with the quark-hadron duality phenomenon to occur, factorization may still be possible [tj [ToL fl~fl| . 

The experiment (E00-108) ran in the summer of 2003 in Hall C at Jefferson Lab. An electron beam with a 
current ranging between 20 and 60 fiA was provided by the CEBAF accelerator with a beam energy of 5.5 GeV. 
Incident electrons were scattered from a 4-cm-long liquid hydrogen or deuterium target and detected in the Short 
Orbit Spectrometer (SOS). The SOS central momentum remained constant throughout the experiment, with a value 
of 1.7 GeV. The electroproduced mesons (predominantly pions) were detected in the High Momentum Spectrometer 
(HMS), with momenta ranging from 1.3 to 4.1 GeV. The experiment consisted of two parts: i) at a fixed electron 
kinematics of (x, Q 2 ) = (0.32, 2.30 GeV 2 ) the central HMS momentum was varied to cover a range of 0.3 < z < 1.0; 
and ii) similarly, at z = 0.55, the electron scattering angle was varied, at constant momentum transfer angle, to span 
a range in x from 0.22 to 0.58. Note that this corresponds to an increase in Q 2 , from 1.5 to 4.2 GeV 2 . The invariant 
mass squared, W 2 , is typically 5.7 GeV 2 and always larger than 4.2 GeV 2 , well in the deep inelastic region, and all 
measurements were performed for both ir + and tt~ . 

Events from the aluminum walls of the cryogenic target cell were subtracted by performing substitute empty target 
runs. Scattered electrons were selected by the use of both a gas Cherenkov counter and an electromagnetic calorimeter. 
Pions were selected using the coincidence time difference between scattered electrons and secondary hadrons. In 
addition, an aerogel detector was used for pion identification [l2T |. For kinematics with pion momenta above 2.4 GeV 
a correction was made to remove kaons from the pion sample, 10% in the worst case (at z ~ 1), as determined from 
the electron-hadron coincidence time. From a measurement detecting positrons in SOS in coincidence with pions in 
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HMS, we found the background originating from ir° production and its subsequent decay into two photons and then 
electron-positron pairs, negligible. 

We modelled semi-inclusive pion electroproduction [13T |, following the high-energy expectation of Eq. @. We 
used the CTEQ5 next-to-leading-order (NLO) parton distribution functions to parameterize q(x,Q 2 ) [lj], and the 
fragmentation function parameterization for (z, Q 2 ) + D' q ^(z,Q 2 ), with D + {D~) the favored (unfavored) 

fragmentation function, from Binnewies et al. [151 ]. The remaining unknowns are the ratio of D~ / D + , taken from a 
HERMES analysis fig , the slope b of the px dependence, and the factors A and B describing the (j) dependence. 

We can not constrain b well within our own data set due to the limited {pt,4>) acceptance of a magnetic spectrometer 



setup. Here, with the possible strong correlation between the pr and (j) dependence [17J], additional assumptions are 
required. Hence, we will use the slope b from an empirical fit to the HERMES pr dependence (b i=s 4.66 GeV -2 ) 
[18j |. Our own best estimate is b — 4.0 ± 0.4 GeV -2 , with no noticable differences between 6- values extracted from 
the pT-dependence of either ir + and ir~ data, or 1 H and 2 H data, somewhat lower than the HERMES slope. We do 
find a (j) dependence in our data, with typical parameters of A = 0.16 ± 0.04, and B — 0.02 ± 0.02, for an average 
< pt > = 0.1 GeV. These 0-dependences become smaller to negligible in the ratios of cross sections shown later. 
Similarly, we find a Q 2 -dependence in our data that differs from the factorized high-energy expectation, but this 
does not affect the results shown below. Of course, these findings do cast doubt on the strict applicability of the 
high-energy approximation for our experiment. 

Within our Monte Carlo package, we estimated two non-trivial corrections to the data. Radiative corrections were 
applied in two steps. We directly estimated the radiation tails within our semi-inclusive pion electroproduction data 
using the Monte Carlo. In addition, we explicitly subtracted radiation tails coming from the exclusive reactions 
e + p — ► e' + 7r + + n and e + n — > e' + ir~ + p. For these processes, we interpolated between the low-IF 2 , low-Q 2 



predictions using the MAID model [19] and the higher- W data of Brauel et al. and Bebek et al. [2fJ, [2l|. We 



subtracted events from diffractive p production, using PYTHIA [22j to estimate the p(e,e'/9°)p cross section with 



similar modifications as implemented by the HERMES collaboration [18!, l23| . We also made a 2% correction to the 
deuterium data to account for the loss of pions traversing the deuterium nucleus [24| . 

The 1,2 H(e,e / 7r ± )A cross sections as measured at x = 0.32 are compared with the results of the simulation in Fig.[TJ 
as a function of z. The general agreement between data and Monte Carlo is excellent for z < 0.65. Within our 
kinematics (pr ~ 0), M 2 is almost directly related to z, as M 2 w M 2 + Q 2 (l/x — 1)(1 — z). Hence, the large excess 
at z > 0.8 in the data with respect to the simulation mainly reflects the 7 TV — 7rA(1232) transition region. Indeed, in 
e.g. a typical 1 H(e,e'7r _ )AT spectrum one can see one prominent A(1232) resonance, and only some small structure 
beyond [13, [HJ. Apparently, above M 2 — 2.5 GeV 2 or so, there are already sufficient resonances to render a spectrum 
mimicking the smooth ^-dependence as expected from the Monte Carlo simulation following the factorization ansatz 
of Eq. ^j. Lastly, the fast drop of the simulations at large z may be artificial. Whereas fragmentation functions have 
been well mapped up to z — 0.9 at the LEP collider (2||, to better than 50%, there remain questions for semi-inclusive 
pion production at lower Q 2 . Here, the fragmentation functions could well flatten out [261 ] . as also included in the 
Field and Feynman expectations [271 ]. that tend to produce more particles at lower energies beyond z = 0.7 or so. 

To quantify the surprising resemblance of semi-inclusive pion electroproduction data in the nucleon resonance region 
with the high energy prediction of Eq. @, we formed simple ratios of the measured cross sections, insensitive to the 
fragmentation process (assuming charge symmetry) at leading order (LO) in a s . If one neglects strange quarks and 
any pr-dependence to the parton distribution functions, these ratios can be expressed in terms of u and d parton 
distributions, as follows 

<Jp(ir + ) + er p (7r~) _ 4:u(x) + 4u(x) + d(x) + d(x) 



Vd{K + ) + Odi^ ) 5(u(x) + d{x) + u{x) + d(x)) ' 
Q'p(7r + ) - Vp{K~) _ ^u v (x) - d v (x) 
Vd{~K+) - a d (-K-) 3(u v (x) + d v (x)) 



(4) 



with the notation a p (ir + ) refering to the ir + pion electroproduction cross section off the proton, u = u v + S, d = d v +d, 
and the Q 2 -dependence left out of these formulas for convenience. These ratios allow us to study the factorization 
ansatz in more detail, with both ratios rendering results independent of z (and pt)- 

We show our results in Fig. [2j with the solid (open) symbols reflecting the data after (before) subtraction of the 
diffractive p contributions. The hatched areas in the bottom indicate the estimated systematic uncertainty. The 
shaded bands reflect the expectations under the assumptions described above (factorization, no strange quark effects, 
charge symmetry for the fragmentation functions), and include a variety of calculations, using both LO and NLO 
(MS and valence) parton distribution functions from the GRV collaboration, and NLO calculations from the CTEQ 
collaboration 
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Our data are remarkably close to the near-independence of z as expected in the high-energy limit, with the clearest 
deviations in the region of z > 0.7, approaching on the A(1232) residual mass region. Within 10% we find perfect 
agreement beyond this region. 

Using the deuterium data only, the ratio of unfavored to favored fragmentation functions D~ / D + can be extracted. 
This ratio is, to a good approximation, at LO simply given by 

(4-^)/(4^-l). (5) 

In the high-energy limit, this ratio should solely depend on z (and Q 2 ), but not on x. The results are shown in Fig.[3j 
with the closed (open) symbols reflecting the data after (before) subtraction of the diffractive p contributions. The 
solid curves are a fit to the HERMES data for the same ratio The dashed curve is the expectation (1 — z)/(l + z) 
according to Field and Feynman for independent fragmentation [27|]. The hatched areas indicate the systematic 
uncertainties, dominated by uncertainties due to the two non-trivial corrections discussed above. 

We observe that the extracted values for D~/D + closely resemble those of the HERMES experiment [la ]. The 
data show a near-independence as a function of x, as expected from Eq. ((21), and a smooth slope as a function of z, 
reflecting a fit to the higher-energy HERMES data, all at > 4 GeV 2 . This is quite remarkable given that our data 
cover the full resonance region for the residual system X, M 2 < M 2 < 4.2 GeV 2 . Apparently, there is a mechanism 
at work that removes the resonance excitations in the ir + /ir~ ratio, and hence the D~ / D + ratio. We note that both 
our data and the fit to the higher-energy HERMES data far exceed the Field and Feynman expectations for large z. 

The mechanism above can be simply understood in the SU(6) symmetric quark model. Close & Isgur [Tl| applied 
this to calculate production rates in various channels in semi- inclusive pion photoproduction, 7 V — > ttX. The pattern 
of constructive and destructive interference, which was a crucial feature of the appearance of duality in inclusive 
structure functions, is in this model also repeated in the semi-inclusive case. The results suggest an explanation 
for the smooth behavior of D~ /D + = /D* for a deuterium target in Fig. [3] The relative weights of the 
photoproduction matrix elements, summed over p and n, is for tt + production always 4 times larger than for 7r~ 
production. In the SU(6) limit, therefore, the resonance contributions to the ratio of Eq. ([5|) cancel exactly, leaving 
behind only the smooth background, as would be expected at high energies. This may account for the glaring lack of 
resonance structure in the resonance region fragmentation functions in Fig. [3l 

In summary, we have measured charged-pion (ir^) electroproduction cross sections for both hydrogen and deuterium 
targets. Our data cover the region where the missing mass of the residual system X is in the resonance region. We 
observe for the first time the quark-hadron duality phenomenon in such reactions, in that such data equate the high- 
energy expectations. We have quantified this behavior by constructing several ratios from these data, that exhibit, at 
low energies, the features of factorization in an electron-quark scattering and a subsequent quark-pion fragmentation 
process. Furthermore, the ratio of favored to unfavored fragmentation functions closely resembles that of high energy 
reactions, over the full range of missing mass. This observation can be explained in the SU(6) symmetric quark 
model. 
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± 





.0322 





.5131 


± 





.0180 


1 


.4767 


± 





.0459 





.9665 


± 





.0979 


0. 


.685 


0. 


.9577 


± 





.0320 


0. 


.5469 


± 





.0191 


1 


.4249 


± 





.0432 





.9206 


± 





.0923 


0. 


.699 


0. 


8742 


± 





.0296 





.4948 


± 





.0177 


1 


.3854 


± 





.0427 





.9392 


± 





.0925 


0. 


.713 


0. 


.8428 


± 





.0286 





.4895 


± 





.0177 


1 


.3738 


± 





.0423 





.9225 


± 





.0924 


0. 


.727 


0. 


.8349 


± 





.0283 


0. 


.4788 


± 





.0171 


1 


.4706 


± 





.0431 





.8840 


± 





.0887 


0. 


741 


0. 


7730 


± 





.0261 





.4720 


± 





.0165 


1 


.3586 


± 





.0396 





.8252 


± 





.0821 


0. 


.755 


0. 


,7029 


± 





.0222 


0. 


.4195 


± 





.0147 


1 


.2689 


± 





.0345 





.7955 


± 





.0732 





.769 


0. 


,6755 


± 





.0216 


0. 


.4409 


± 





.0156 


1 


.2450 


± 





.0336 





.7247 


± 





.0674 


0. 


783 


0. 


,5922 


± 





.0192 





4336 


± 





.0145 


1 


.2769 


± 





.0340 





.7418 


± 





.0674 


0. 


.797 


0. 


,5987 


± 





.0194 





.4654 


± 





.0155 


1 


.3623 


± 





.0360 





.7660 


± 





.0687 


0. 


.811 


0. 


,6461 


± 





.0192 


0. 


.5136 


± 





.0163 


1 


.4704 


± 





.0375 





.7862 


± 





.0692 


0. 


.825 


0. 


6808 


± 





.0199 


0. 


.5425 


± 





.0166 


1 


.5252 


± 





.0375 





.7965 


± 





.0697 





.839 





7641 


± 





.0213 


0. 


.5979 


± 





.0177 


1 


.6375 


± 





.0397 





.8278 


± 





.0698 


0. 


.853 


0. 


,7667 


± 





.0212 


0. 


.6530 


± 





.0186 


1 


.6277 


± 





.0369 





.8383 


± 





.0670 


0. 


.867 


0. 


,8048 


± 





.0210 


0. 


.6553 


± 





.0171 


1 


.6313 


± 





.0341 





.7840 


± 





.0598 





.881 





,6894 


± 





.0212 


0. 


.4966 


± 





.0144 


1 


.4526 


± 





.0311 





.7161 


± 





.0574 


0. 


.895 


0. 


,5098 


± 





.0209 


0. 


.3664 


± 





.0131 


1 


.0932 


± 





.0279 





.4969 


± 





.0500 


0. 


.909 


0. 


,2692 


± 





.0232 


0. 


.2400 


± 





.0122 





.7644 


± 





.0274 





.3685 


± 





.0512 


0. 


.923 





0222 


± 





.0350 


0. 


.1380 


± 





.0120 





.5219 


± 





.0365 





.2641 


± 





.0630 



TABLE I: Cross sections as a function of z in nb/GeV 3 /sr coresponding to Fig. 1. Errors are statistical only. 
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FIG. 1: The 1 ' 2 H(e,e , 7r ± )X cross sections at x = 0.32 as a function of z in comparison with Monte Carlo simulations (dashed 
curves) starting from a fragmentation ansatz (see text). The various cross sections have been multiplied as indicated for the 
purpose of plotting. See Table I for numerical values. 



8 



+ 



to 

+ 



to' 

to 
i 



1.6 



14 



3 12 



0.8 

3 
2.5 
2 

1.5 
l 

0.5 





'irfffftttif/J /£///////////, 



-}} " H 



^ +* \k A 



0.4 



0.5 



0.6 



0.7 



0.8 



FIG. 2: The ratio of proton to deuterium results of the sum (top) and difference (bottom) of ir + and 7r~ cross sections as a 
function of z, at a; = 0.32. Closed (open) symbols reflect data after (before) events from coherent p production are subtracted 
(see text). The symbols have been sightly offset in z for clarity. The hatched areas in the bottom indicate the systematic 
uncertainties, whereas the shaded bands represent a variety of calculations, at both leading order and next-to-leading-order of 
a s , of the shown ratio OH^. See Tables II and III for numerical values. 
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FIG. 3: Top: The ratio of unfavored to favored fragmentation function D~ /D + as a function of x at z = 0.55, evaluated at 
leading order of a 3 from the deuterium data. The triangles (square) reflect HERMES (EMC) data [Isl . [29| . with the solid 
curve a fit to HERMES data. Further symbols and the hatched area are as in Fig. 2. Bottom: Same as top, but now as a 
function of z for x — 0.32. The dashed curve represents the expectation [27| under the independent fragmentation hypothesis. 
See Tables IV and V for numerical values. 
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z R+ (after p) R+ (before p) 






.349 


1 


.1790 


± 





,0750 


1. 


.1670 


± 





.0710 





363 


1 


.2700 


± 





.0700 


1 


,2570 


± 





.0650 





.377 


1 


.1180 


± 





.0550 


1 


,1060 


± 





.0510 





.391 


1 


.1850 


± 





.0510 


1. 


,1700 


± 





.0480 





.405 


1 


.2110 


± 





.0500 


1 


,1970 


± 





.0460 





.419 


1 


.1040 


± 





.0480 


1. 


,0900 


± 





.0440 





.433 


1 


.2070 


± 





.0460 


1. 


,1910 


± 





.0420 





.447 


1 


.1490 


± 





.0450 


1 


1320 


± 





.0410 





.461 


1 


.2330 


± 





.0490 


1. 


,2120 


± 





.0440 





.475 


1 


.1960 


± 





.0510 


1. 


,1730 


± 





.0450 





.489 


1 


.2090 


± 





.0500 


1 


,1880 


± 





.0440 





.503 


1 


.1340 


± 





.0460 


1 


,1140 


± 





.0400 





.517 


1 


.2660 


± 





.0520 


1. 


,2390 


± 





.0460 





.531 


1 


.2300 


± 





.0530 


1 


,2040 


± 





.0470 





.545 


1 


,2890 


± 





.0540 


1. 


,2600 


± 





.0470 





.559 


1 


.2010 


± 





.0530 


1 


,1720 


± 





.0460 





.573 


1 


.2720 


± 





.0490 


1 


,2390 


± 





.0420 





.587 


1 


.2890 


± 





.0490 


1 


,2540 


± 





.0420 





.601 


1 


.1760 


± 





.0450 


1 


,1430 


± 





.0380 





.615 


1 


.2700 


± 





.0480 


1 


,2340 


± 





.0410 





.629 


1 


.1410 


± 





.0430 


1 


,1110 


± 





.0370 





.643 


1 


.1370 


± 





.0420 


1 


,1070 


± 





.0360 





.657 


1 


.1760 


± 





.0400 


1. 


.1450 


± 





.0340 





.671 


1 


.1690 


± 





.0420 


1 


.1370 


± 





.0350 





.685 


1 


.1870 


± 





.0410 


1 


,1570 


± 





.0350 





.699 


1 


.1450 


± 





.0400 


1 


,1140 


± 





.0340 


n 

\J 


71 ^ 


I 


1 TOO 




n 

\J 


run n 


I 


i ^nn 




n 


■ UOOU 





.727 


1 


.1050 


± 





.0390 


1. 


,0800 


± 





.0330 





.741 


1 


.1220 


± 





.0380 


1. 


.0980 


± 





.0330 





.755 


1 


.0580 


± 





.0340 


1 


0360 


± 





.0300 





.769 


1 


.0650 


± 





.0340 


1 


,0480 


± 





.0300 





.783 





,9990 


± 





.0320 





,9870 


± 





.0280 



TABLE II: The ratio R^ d of proton to deuterium results of the sum of ix + and it cross sections as a function of z, at 
x = 0.32. Results reflect data after (left) or before (right) events from coherent p production are subtracted (see text). 
Errors are statistical only. The systematic error is given by 0.0854 - 0.5462 + I.I82 2 - 0.778z 3 for 0.3 < 2 < 0.68 and by 
-3.313 + 13.432 - 17.99z 2 + 7.995z 3 for 0.68 < z < 0.80. 
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z R pd (after p) R pd (before p) 






.363 


1 


.4440 


± 


0. 


.2710 


1 


.4570 


± 


0. 


2730 





.405 


2 


.0860 


± 


0. 


3190 


2 


.0910 


± 


0. 


3210 





.447 


1 


.6750 


± 


0. 


.2180 


1 


.6300 


± 


0. 


2120 





.489 


1 


.5690 


± 





2080 


1 


.5210 


± 


0. 


.2030 





.531 


1 


.8670 


± 


0. 


2290 


1 


.8610 


± 


0. 


.2280 


o 


.573 


I 


.9340 


± 


0. 


2220 


x 


.9720 


± 


0. 


.2280 





.615 


1 


.7160 


± 


0. 


1770 


1 


.7640 


± 


0. 


1830 





.657 


1 


.4470 


± 


0. 


1440 


1 


.4930 


± 


0. 


1500 





.699 


1 


.4260 


± 


0. 


.1540 


1 


.5070 


± 


0. 


1680 





.741 


1 


.0530 


± 


0. 


.1120 


1 


.0740 


± 


0. 


1180 





.783 





.5590 


± 


0. 


.0780 





.5270 


± 


0. 


0800 



TABLE III: The ratio R~ d of proton to deuterium results of the difference of tt + and 7r~ cross sections as a function of z, at x 
= 0.32. Results reflect data after (left) or before (right) events from coherent p production are subtracted (see text). Errors 
are statistical only. The systematic error is given by 0.046 — 0.203z + 0.538.Z 2 — 0.163z 3 . 



z D-/D+ (after p) D-/D+ (before p) 






.342 





,4620 


± 


0. 


.0710 





.4920 


± 


0. 


0700 





.370 





4196 


± 


0. 


0475 





.4449 


± 





.0465 





.398 


0. 


.4838 


± 


0. 


.0453 





.5126 


± 


0. 


.0438 





.426 





4764 


± 


0. 


.0429 





.5087 


± 


0. 


0411 





.454 





4575 


± 





0414 





.4940 


± 





0392 





.482 


0. 


.4425 


± 





.0413 





.4837 


± 


0. 


.0395 





.510 





,4059 


± 


0. 


.0318 





.4530 


± 


0. 


0306 





.538 





.3635 


± 


0. 


.0270 





.4134 


± 





0257 





.566 





3699 


± 





0266 





.4288 


± 





0253 





.594 


0. 


,3638 


± 


0. 


.0274 





.4280 


± 


0. 


,0267 





.622 





,3448 


± 


0. 


.0298 





.4124 


± 





,0284 





.650 





3157 


± 





.0289 





.3853 


± 





,0279 





.678 


0. 


3587 


± 


0. 


.0314 





.4376 


± 





,0307 





.706 


0. 


,3934 


± 


0. 


.0327 





.4800 


± 





0319 





.734 





3137 


± 





0273 





.3889 


± 





,0264 





.762 





3164 


± 





.0254 





.3911 


± 





,0254 





.790 





,2738 


± 


0. 


.0223 





.3375 


± 





,0223 





.818 





,2625 


± 





0198 





.3177 


± 


0. 


,0198 





.846 





,2380 


± 





.0177 





.2808 


± 


0. 


,0168 





.874 





,2294 


± 


0. 


.0161 





.2607 


± 





0159 





.902 





,2423 


± 





0243 





.2555 


± 





,0238 





.930 





,3025 


± 





.0739 





.2944 


± 


0. 


,0724 



TABLE IV: The ratio of unfavored to favored fragmentation function D~ / D + as a function of z at x = 0.32, evaluated at 
leading order of a s from the deuterium data. Results reflect data after (left) or before (right) events from coherent p production 
are subtracted (see text). Errors are statistical only. The systematic error is given by 2.445 - 18.434,2 + 46.196z 2 - 38.194z 3 
for 0.2 < 2 < 0.3, 0.273 - 1.684z + 3.461z 2 - 2.132z 3 for 0.4 < z < 0.835, and -24.15 + 89.532.3 - 110.09z 2 + 44.974z 3 for 
0.835 < z < 0.935. 
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x D~/D+ (after p) D~/D+ (before p) 



n 

u 


.213 


n 

u 






n 


UOOJ 


n 
u 


^89 


j- 


n 
u 


nsoo 

uouu 


n 
u 




n 
u 


4979 


_l_ 


n 
u. 




n 
u 


478Q 


ZZ\ 


n 
u. 




n 
u 




n 
u 


4008 

.^lUUO 


m 


n 

u. 


n^?4i 


n 

u 


4479 


-L 

zn 


u 




n 

u 


987 

.ZO 1 


n 
u 




-J- 


n 


n^i i 


n 
u 




-J- 


n 
\j 


D9Q4 


n 
u 






4D4Q 




n 


D98Q 
Uioy 


n 
u 


444R 


-L 
zz\ 


n 
u. 


D977 
uz / / 


n 

u 


.OOO 


n 

u 


4978 

.^±Z i o 


-1- 


n 


098^ 

UZOJ 


n 
u 


1UUU 


ZZ 


n 

u 


0974 


n 


-OUO 


n 


.0004: 




n 


0252 


n 
u 


-OUO 1 




n 
u. 


0242 


n 


-OOO 


n 




-1- 


n 




n 
u 


.oyo i 


-1- 


n 


09^ 


n 


4.1 ^ 


n 

u 


^47fi 


-1- 


n 


0262 


n 
u 


^7^9 


-1- 


n 


094Q 


n 


4^8 


n 

u 


^!Q14 


zz 


n 




n 
u 


41 77 


-1- 


n 
u 


0987 





.463 





.3907 


± 





.0320 





.4142 


± 


0. 


.0310 





.488 





.4198 


± 


0. 


.0362 





.4420 


± 





0349 





.513 





.4336 


± 





0403 





.4546 


± 





.0395 





.538 





.4202 


± 





.0454 





.4385 


± 


0. 


.0440 





.562 





.4721 


± 


0. 


.0581 





.4890 


± 





0572 





.588 





.3533 


± 





.0553 





.3668 


± 


0. 


.0539 



TABLE V: The ratio of unfavored to favored fragmentation function D~ /D + as a function of x at z = 0.55, evaluated at leading 
order of a s from the deuterium data. Results reflect data after (left) or before (right) events from coherent p production are 
subtracted (see text). Errors are statistical only. The systematic error is given by 0.126 - 0.669x + 1.445x 2 - 1.119a; 3 . 



